Solvation-induced one-dimensional polarons and electron transfer 
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When a one- dimensional (ID) semiconductor nanostructure is immersed in a sluggish polar sol- 
vent, fluctuations of the medium may result in the appearance of localized electronic levels inside 
the band gap. An excess charge carrier can occupy such a level and undergo self-localization into a 
large-radius adiabatic polaron surrounded by a self-consistent medium polarization pattern. Within 
an appropriately adapted framework of the Marcus theory, we explore the description and qualita- 
tive picture of thermally activated electron transfer involving solvation-induced polaronic-like states 
by considering transfer between small and ID species as well as between two ID species. Illustrative 
calculations are performed for tubular geometries with possible applications to carbon nanotube 
systems. 

PACS numbers: 31.70.Dk, 71.38.-k, 82.20.Yn, 82.45.Yz 



I. INTRODUCTION 



Electron transfer reactions mediated by fluctuations 
of the (polar) environment are of great importance for 
a variety of physico-chemical and biological processes. 
Descriptions of the Marcus theory of electron transfer 
and various elaborations on the topic are staples in nu- 
merous textbooks (e. g, 0, 0, S i, ij). In terms of 
traditional application areas, one may distinguish elec- 
tron transfer between small species (chemical reactions in 
solutions), as schematically illustrated in Fig. [lja), and 
transfer between small species and bulk electrodes (elec- 
trochemistry), Fig. [He). In the former case, the trans- 
fer occurs between two discrete electronic states; in the 
latter case the transfer is between a discrete level and a 
rigid continuous band of electronic states. The subject of 
our interest in this paper is thermally activated electron 
transfer to/from one-dimensional (ID) semiconducting 
nanostructures, as sketched in Fig. [IJb) and exemplified 
by systems like conjugated polymers [g], nanotubes and 
nanowires 0- ID systems are known to be particularly 
sensitive to potential energy fluctuations resulting in the 
formation of bound states below the continuum even in 
shallow potential wells |8j. In fact, when an excess charge 
carrier is added to a ID semiconductor, the long-range 
Coulomb interaction with the surrounding polar solvent 
may result in self-trapping of the carrier, where a local- 
ized electronic state is accompanied by a self-consistent 
dielectric polarization pattern of the medium. Forma- 
tion of such solvation-induced ID polarons has been dis- 
cussed recentl y b oth within simplified theoretical models 
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1 12l.Tl3|] and at the level of ab initio computa- 
Ha ]. The sketch of electron transfer in Fig.QJb) 
involving polaron-like localized electronic states is thus, 
in a sense, intermediate between traditional pictures of 
Fig. (Ha) and (c). 

It should be noted here that ID semiconductor nanos- 
tructures in contact with polar solvents are candidates for 
various technological applications involving fundamental 
redox reactions, such as in photoelectrochemistry, energy 
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FIG. 1: Schematic pictures of electron transfer (a) between 
two small species, (b) between a small and a ID species, and 
(c) between a small species and a bulk semiconductor. The 
second row depicts electronic levels (and bands) and their 
modulation by medium fluctuations to promote iso-energetic 
electron transfer. 

harvesting and sensing (e. g., @, d [H El El HI HI, 
[22|). As we will be using later specific illustrative cal- 
culations for small-diameter tubular structures, semicon- 
ducting single-wall carbon nanotubes (SWCNTs) are a 
particularly widely known example of such. Redox chem- 
istry of CNTs has been deemed an "emerging field of 
nanoscience" [23[ and solvato chromic effects in CNTs are 
being intensely researched [HI • 

Considerations of the electron transfer problem [H,[l,[l, 
Hill involve various aspects and regimes. In this paper we 
are concerned only with a very limited scope of questions. 
The regime we are interested here is usually referred to as 
high-temperature (T) , non-adiabatic transfer epitomized 
by the famed Marcus expression (e.g., [1,0 ) 
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for the transfer rate between discrete electronic states 1 
("donor") and 2 ("acceptor"). Here t2\ is the electronic 
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transition matrix element, while the activation energy 
(£ r2 — £ t \ + E r ) 2 



(2) 



is determined by the difference of the fully equilibrated 
(free) energies £ t2 and £ t \ of the respective states, and 
the so-called reorganization energy E r . 

We restrict our consideration to the outer sphere reor- 
ganization referring to the interaction of the charge car- 
rier with the polarization P(r) of the polar medium. (For 
a related discussion of the effects of the interaction with 
ID lattice distortions in the context of the electron in- 
jection into a deformable polymer chain, see Ref. [HI].) 
Expression J]} can be derived from the averaging of the 
Fermi golden rule 
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(3) 

for the electronic transition over various spatial distri- 
butions {S>~P) of harmonic solvent fluctuations under as- 
sumptions of the Boltzmann probability ^i(P) with the 
electron in state 1, linear coupling of purely electronic 
energies E\ and E 2 to fluctuations, and t 2 i being inde- 
pendent of P. Equation ([3]) very nicely reflects the spirit 
of the Marcus theory, where iso-energetic electron trans- 
fer takes place in configurations created by appropriate 
solvent fluctuations. In what follows we discuss applica- 
tions of this approach to electron transfer between a small 
species and a ID semiconductor as well as between two 
parallel ID structures. Our emphasis would be on the 
qualitative picture of thermally activated transfer and on 
the evaluation of the corresponding activation energies. 

While, for certainty, we explicitly discuss the transfer 
of an electron in the context of the conduction band of a 
ID semiconductor, it should be absolutely clear that the 
same picture applies to the transfer of a hole in the con- 
text of the valence band of the semiconductor. One could 
as well use the language of the electronic excitation ener- 
gies equally applicable to both electrons and holes. Also, 
the principle of detailed balance in thermal equilibrium 
allows one to immediately relate rates ki^ 2 of the for- 
ward ("to") and k 2 ^i of the reverse ("from") processes 
0, Q and their respective activation energies. 



II. ELECTRON SOLVATION ON SMALL 
SPECIES AND ON ID SEMICONDUCTOR 

As an introduction to the further discussion, in this 
section we compare solvation of an electron (charge q) 
separately on a small spherical species and on a ID semi- 
conductor structure. The question of interest here is a 
relationship between the (lowest) purely electronic en- 
ergy E as affected by the solvent polarization P(r) and 
the minimal possible total free energy £ of the electron- 
solvent system for a given E. The total energy includes, 
besides E, the (free) energy U stored in the solvent po- 
larization. Our question of interest can also be expressed 



as that of minimal possible energy U for a given E. The 
bare, unaffected by the polarization, electronic energy 
will be denoted as E°. 

We find it very convenient to represent the state of 
solvent polarization P(r) by a distribution of fictitious 
charges Zi that would cause this polarization pattern as 
an equilibrium dielectric response 0, licj . In terms of 
such charges, the energy stored in the polarization would 
then be generically given as 
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(4) 



while the contribution to the electronic potential energy 
described by 



(5) 



The interaction kernel elements Vij and V e i are deter- 
mined by solutions of the corresponding electrostatic po- 
tential problem for specific spatial distributions (and po- 
sitions) of fictitious charges Zi and the real charge q. 

One recalls that polarization P relevant for the elec- 
tron transfer problem is not the total solvent polarization 
but its slower (orientational) component [H, HI HI EL HI- 
For a polar environment, such as the Debye solvent, char- 
acterized by two dielectric constants: static e s and high- 
frequency eoo, the effective dielectric constant e* is usu- 
ally defined via 

1/e* = l/ £oo - l/e s 

to represent the dielectric response due to the slower 
polarization component. (For typical [J, [l(| solvents, 
e s ^ £00 so that e* ~ €00.) The fast component of the 
total polarization acts "instantaneously" in the present 
context, and its effect is considered included in the defi- 
nition of the bare electronic energies E°. For simplicity, 
we restrict our attention to the uniform dielectric envi- 
ronment, see Refs. [13; f° r a discussion of the effects 
of various dielectric conditions. 

The well-known case of a small spherical species of ra- 
dius r is very simple as it corresponds to the fixed on- 
the-sphere distributions for both charge q and single fic- 
titious charge z. That results in Vn = V e \ = l/e*r. The 
electronic energy in this case is given entirely by Eq. ([5]): 



E{z} = E° -qz/e* 
while the polarization energy ((4]) by 
U{z] = z 2 /2e*r. 



(6) 



(7) 



The total system energy £{z} = E{z}+U{z} is a familiar 
parabola achieving its minimum of £ x = E — E T at z = q 
with the individual reorganization energy 



E T = q 2 /2e*r. 



(8) 
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FIG. 2: Scaled functional dependences of the optimal (a) po- 
larization energy U and (b) total energy £ for a given purely 
electronic energy E. Equlibrated solvation corresponds to the 
minima of total energies in panel (b). Long-dash lines show 
the benchmark parabolic behavior for small species; solid lines 
describe solvation on a ID species with the confinement pa- 
rameter clb/R = 4, short dashes are for the increased confine- 
ment of as/R — 10. 



As is evident from Eqs. ([6])- (JSj , the question we posed in 
the beginning of this section is answered by the parabolic 
dependence 



U(E) = {E- E y/4E I 



(9) 



for the small-species case. These benchmark dependences 
of U{E) and £{E) = E + U(E) are shown in Fig. [2] in 
terms of the energy scale © • 

Description of the adiabatic electron solvation on a ID 
semiconductor is more involved as both the real (elec- 
tron) charge q and the fictitious charge z can have vari- 
able spatial extents along it (coordinate x): 



q = dxp(x), z — dxp*{x), 



(10) 



where p(x) and p*(x) are the respective linear charge 
densities. The former is determined by the (normalized) 
electron wave function ip(x): 



following from the solution of the Schrodinger equation 
H(x)i/}(x) = (E - E°)ip(x). (12) 
The Hamiltonian 



h 2 d 2 
2m ox z 



(13) 



here is composed of the kinetic energy with the effec- 
tive mass m and the potential energy due to the solvent 
polarization: 



dy V(x - y) p f (y). 



(14) 



This potential energy (|14[) is of course just a continuum 
version of Eq. ([5]) , while the corresponding energy of the 
polarization ([¥]) is now given by 



U{p f } 



dx dy pf (x) V (x — y) pf (y) . (15) 



An explicit example of the kernel V(x) is discussed below. 

It is understood that only the lowest energy eigenvalue 
E < E° is of relevance for us here. We now want to 
find configurations with the minimal possible energy U 
([TB|l that would yield a given value of E. The constrains 
imposed are conveniently taken into account in the min- 
imization of the auxiliary functional 

U = U + a ((H) - (E - E )) + fitf, (16) 
where a and (3 are the Lagrange multipliers, 



(H) 



dx tj)* (x)H(x)?p(x) 



and 



N = 



J dx\^(x)f 



(17) 



One immediately finds that the variation of U in Eq. (jT6|) 
over ip*(x) would lead to Eq. (T2J) for (3 = a(E° - E). 
The variation over p* (x), on the other hand, establishes 
that 



(18) 



p(x) = q\ip(x)\ 2 



(11) 



p f {x) = ap(x), 



with p(x) defined by Eq. (fTTj) . Comparison with Eq. (JTT 
shows that multiplier a has the meaning of the relative 
magnitude of the fictitious charge: a — z/q. 

Using a parametrically, one can now solve a self- 
consistent problem of Eqs. ![T2"]). (fTTj) , (fT5|). and (H3J) iter- 
atively to find the required dependences of U{E) and 
£{E) = E + U(E). The minimum of E = E° - E r 
is achieved at a = 1 - this is the equilibrium polaron 
state. The individual reorganization energy E r in this 
case is called the polaron binding energy, J3&. We have 
already discussed at length [l(| [H[ the spatial extent and 
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the binding energy of equilibrium solvation-induced po- 
larons; for a practical range of parameters, the scales 
of those quantities are respectively established by the 
Bohr radius and Rydberg energy for the corresponding 
3D Coulomb problem: 



as = e*h 2 /mq 2 , Ry = q 2 /2e*aB- 



(19) 



For a perspective, if one were to take representative val- 
ues of e* = 3 and of the effective mass of 0.05 m e (that 
would correspond to SWCNT of radius R ~ 8 A [2$), 
then Eq. (19} would result in a B ^ 32 A and Ry ~ 76 
me V. A twice larg er effective mass (say, in a smaller ra- 
dius SWCNT [26j) would decrease as and increase Ry 
by a factor of 2; large variations can also be induced by 
changes in e*. 

Figure [2] shows an example of numerical results for 
U(E) and £{E) scaled with respect to polaron's E T = 
Eb- The deviations from the parabolic dependences © 
of the small species are evident. As it frequentl y t akes 
place for excitonic and polaronic problems in ID jl2l . [l3l . 
I271 . [28| , a power-law turns out to be a good approximate 
representation in this case as well: 

U/Er ^ 0.3((E° - E)/E t ) 1 - 7 . 

This approximate relationship is, of course, not universal 
and small parameter-dependent deviations are illustrated 
in the figure. It should also be noticed that the standard 
relationship E T = (E° — E)/2 = U for equilibrated sol- 
vation on small species does not hold for a ID polaron 
(compare minima of the curves in Fig. [2^b)). This is a 
consequence of the additional degrees of freedom (vari- 
able spatial extent) in the polaron case. 

For the numerical calculations above and in the ex- 
amples that follow, we use electrostatic kernels corre- 
sponding to tubular (nanotube) charge distributions of 
transverse radius R. As is well-known [29|], continuum 
models with the purely ID Coulomb interaction lead to 
diverging excitonic and polaronic bindings. Taking into 
account the transverse distribution of the charge (rings 
in the case of the tubular geometry) regularizes that 
behavior, and the binding becomes a gro wing function 
of the confinement parameter as/R (e.g., [U EE H3 ) • 
For nanotubes, the kernel V(x) in Eqs. (fT4"|) and (TT5l) is 
conveniently expressed via its Fourier-transform [IOL |30I | 
V(k) = J dx exp(-ikx) V(x): 



V(k)=2I (kR)K (kR)/e*, 



(20) 



where / and K are the modified Bessel functions appear- 
ing in electrostatic problems with cylindrical symmetry 
31 1 . In earlier work [ll[ we established that the ID po- 
laron picture is valid for as/R > 1; for illustrative cal- 
culations in this paper we use some representative values 
[l3| of the confinement parameter. For specific applica- 
tions to semiconducting SWCNTs, one should note that 
the effective mass in that case scales with tube's radius, 
approximately as m oc 1/R [26], which makes the con- 
finement parameter nearly R- independent. 



III. TRANSFER BETWEEN SMALL SPECIES 
AND ID POLARONIC STATES 

In order to illuminate the commonalities and differ- 
ences with the traditional Marcus setup, we start here 
with a brief exposition of the well-known case of transfer 
between two small spherical species. Our discussion is 
similar to the line of presentation used in Ref. [5j . 

Denoting fictitious charges on donor and acceptor 
species as z\ and z 2 , respectively, the modulation of elec- 
tronic energies in states 1 and 2 in this case is as given 
by Eq. ©: 



Ei — E 1 



Eo — En 



£2 
•'•2 



ri2 



(21) 



(22) 



where r\ and r 2 are the radii of the species and r 12 is the 
distance between their centers. The polarization energy 
flU) is then 



U 



2n 



2r 2 



z\z 2 
ri2 



(23) 



The equilibrated solvation on each of the species corre- 
sponds to the unconstrained minimization of total en- 
ergies £ as discussed in Sec. |TT] with the same results. 
One immediately finds that minimization of £\ = Ei + U 
yields the minimum value of £ r \ = E® — q 2 )2e*r\ at Z\ = 
q, Z2 = 0. The minimization of £2 = E2 + U, on the other 
hand, yields the minimum value of £ r2 = E^ — q 2 /2e*r 2 
at Z\ — 0, z 2 = q. 

While the Marcus result (Eqs. (U) and ©) directly 
follows from employing expressions (p?Tj) - (|2"3"|) in Eq. ([3]), 
here we instead exclusively focus on the electron transi- 
tion point. The optimal state of polarization at the tran- 
sition is again determined by the minimum of the total 
energy but now constrained by the condition E\ = E% 
of equal purely electronic energies (which is the same as 
equal total energies: £± — £ 2 ). Minimization of the aux- 
iliary energy function 



£ = U + E 1 +a(E 2 -E 1 ) 



(24) 



with the Lagrange multiplier a results in relationships 
that will be recurring in what follows: 



Z\ = (1 - a)q, z 2 = aq, Z\ + z 2 



(25) 



The parameter a thus has a meaning of the partitioning 
of the charge q between fictitious charges at the transition 
point. Optimization of Eq. (|24)l yields its value at the 
transition as 



a=(£ l2 -£ I i+E I )/2E I , 
where reorganization energy 

e* \2r\ 2r 2 n 2 



(26) 



(27) 
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consists of donor and acceptor individual reorganization 
energies corrected by the interaction term. As used in 
Eq. ([25]) . parameter a can also play a role of the con- 
venient one-dimensional "reaction coordinate" [5j: from 
a = for the electron in equilibrated state 1 to a = 1 for 
electron in state 2. The resulting parametric dependence 
of total energies in this case is that of famous Marcus 
parabolas: 

£ 1 = £ tl +E I a 2 , £ 2 = £ l2 + E,(l - a) 2 . (28) 

The transition value of a in Eq. ([2"5]) at which £\ = £2 = 
£ tr is of course given by Eq. ([26]) . The activation energy 
of the 1 — > 2 transition: 

E act = £ tI - £ Tl = E r a 2 , (29) 

is precisely the familiar Eq. ([2]). 

The modification of the above approach to the case of 
the acceptor being a ID semiconductor species is done in 
the spirit of the development in Sec. [TIJ where the linear 
distributed charge densities ([TO]) . (fTT|) have been used. 
Particularly, the fictitious charge on the acceptor now 
becomes 




while the energy of the electron on the donor becomes 

Ei=E%- qz l /e*r ~q J dxV 12 (x) p f (x), (30) 

where r is the donor radius, and the interaction kernel 
Vi2 is exemplified below. The analogy between (|30j) and 
earlier used (|2Tj) is evident. 

On the other hand, the energy of the polaronic-like 
electron state on the acceptor, E 2 — E® would now be de- 
termined by the lowest-energy solution of the Schrodinger 
equation (fT2|) with the potential 

y(x) = -qJdy V(x - y) P f (y) - qz 1 V 12 (x). (31) 

We reiterate that the range of parameters we are inter- 
ested in corresponds to the bound electronic states be- 
low the band edge: E 2 < E®. Equation ([31) differs from 
Eq. (Til) by its second term, which represents the inter- 
action with the fictitious charge on the donor (in close 
analogy to a similar interaction term in Eq. {22])). 

Finally, the place of the polarization energy ([23]) in the 
Marcus setup is now taken by 

U = U 1 {z 1 } + U 2 {p f } + z 1 J dxV 12 {x) p f (x) 

where U\{z\} corresponds to the fictitious charge on 
small species ([7]), U 2 {pf} to the linearly distributed fic- 
titious charge (p~5|) . and the last term is their interaction. 
Using, for certainty, a tubular transverse distribution of 
ID charges on a tube of radius R, the intratube interac- 
tions V(x) are described by Eq. (12"0]) while the interaction 



V\ 2 (x) with the small spherical species by the Fourier- 
transform 

V 12 (k) = 2I (kR)K (kd)/e*, (32) 

where d is the distance between the tube axis and the 
center of the sphere. 

Unconstrained (in electronic energy) minimization of 
total energies Ei + U and E 2 +U here results, of course, in 
equilibrated solvation of small species (at z\ = q, z 2 = 0) 
and ID polaron (at z\ — 0, z 2 — q), respectively, as 
described in Sec. [Til The optimal polarization fluctuation 
for iso-energetic electron transfer is constrained by the 
condition E\ = E 2 and found by minimization of the 
auxiliary functional 

£ = U + E X + a(E 2 - Ei) + fW, (33) 

where Lagrange multipliers a and (3 are immediately rec- 
ognized by comparing to Eqs. ([16]) and ([24]) . It then 
comes as no surprise that variations of Eq. (|33[) over 
ip*{x) and pf(x) lead to earlier found results: (3 = 
a(i?2 — E 2 ) for normalization (fl~7]) related coefficient, and 
Eq. (|18p for the relationship between fictitious and real 
charge densities. Furthermore, variation over z\ results 
in z\ = (1 — a)q, and, hence, the same type ([25]) of charge 
q partitioning between fictitious charges is recovered at 
the transition point that was encountered in the conven- 
tional Marcus problem. Evaluating then the difference 
between the total energy £ tr at the optimal fluctuation 
and the equilibrated energy £ T \, one arrives at the acti- 
vation energy for the 1 — > 2 transition as 

£act = + U 2 {p} - q J dxV 12 (x) p(x^j a 2 . (34) 

While the structure of Eq. ([3"4"]) is analogous to the 
structure observed in Eqs. (|29]) . ([27]) . there is a subtle dif- 
ference between them in terms of the dependence on the 
value of the partitioning parameter a at the transition 
point. The parameter E r in Eq. (|29p is independent of a, 
and therefore the a-dependence of the activation energy 
there is strictly parabolic. Differently, the charge density 
p{x) featured in Eq. (|34f) is not the equilibrated polaron 
density (which would be obtained for a = 1) but, rather, 
the self-consistently determined transition state density 
generally dependent on a itself (determined, in turn, by 
the equilibrated energy offset between the species). Thus, 
rigorously speaking, the a-dependence in Eq. ([34]) is not 
strictly parabolic. By the same token, if one were to con- 
sider the 2 — > 1 transition, there would be no strict sym- 
metry between forward and reverse transitions as seen in 
Eq. ([28]) . These subtle features are related to a qualita- 
tive difference in the details of solvation on small and ID 
species (Sec. [TTJ) and reflect the presence of new degrees 
of freedom for the solvated electron in variability of its 
spatial extent on the ID species. 

It is another matter that, quantitatively, the effect due 
to the small species can easily dominate, largely disguis- 
ing those qualitative differences and deviations from the 
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is worthwhile to stress here that the a = point cor- 
responds to the equilibrium polarization pattern of the 
small species (no fictitious charge on the acceptor (|25p). 
Still, the transition occurs into an electronic level below 
the band edge as is clearly seen by comparison with the 
dashed-line rigid-band result. This localized state on the 
acceptor is entirely due to the electrostatic interaction 
between small and ID species as specified by the second 
term in Eq. (|3Tj) . Given this, we cannot exclude that 
the "inverted Marcus regime" [l| could be taking place 
here over some range of energy offsets. To verify this 
conjecture, however, one would have to extend the cur- 
rent framework to explicitly include a multitude of elec- 
tronic states on the ID species, and we do not show it in 
the figure. The reverse 2 — > 1 transition into a discrete 
electronic level of the small species features the inverted 
regime as usual @. 



FIG. 3: Example of calculated activation energies for elec- 
tron transfer between a small (state 1) and a ID (state 2) 
species as a function of the energy offset between the equilib- 
rium energy of the small species and the band edge of the ID 
semiconductor. The dashed lines show the results one would 
obtain for a rigid electronic band unaffected by the solvent 
polarization. See text for details. 



conventional parabolic behavior. Also contributing to 
this "disguise" is a compensatory character of the effects 
related to the second and third terms in parentheses of 
Eq. ([34]) . The former originates in electrostatic interac- 
tions within the ID system, while the latter in interac- 
tions between the small and ID species. 

Figure [3J depicts results of illustrative calculations of 
the activation energies for the case characterized by the 
geometric parameters r = 3.5 A, R — 5 A, d = 10 A, 
and as — 15 A in the corresponding Ry units (fT9|) . 
Shown are activation energies for both forward 1 — > 2 
and reverse 2 — ► 1 transitions. The results calculated 
as described above (solid lines) are compared to what 
one would get if the band states were unaffected by the 
solvent polarization (a rigid band, long-dash lines). In 
the latter case, the reorganization energy E t would be 
entirely due to the small species, and the activation en- 
ergy for the 1 — > 2 transition is well represented by the 
approximation i? act = (E® — £ri + E T ) 2 /AE r (as one can 
numerically confirm from the conventional [f| integration 
of Eq. (p} over band states: 

fcl^band- / dE 2 D{E 2 )k 1 ^ 2 (E 2 ), 

Je° 



with ID density-of-states D{E 2 ) cx l/y/E 2 - E%.) The 
relative displacement of the "center lines" shown in the 
figure indicates the equilibrium polaron binding energy 
Ei, (which is close to 1 Ry in this case). 

Marked in Fig. [3] with "a = 0" and "a = 1" are energy 
offsets where the respective activation energies vanish. It 



IV. TRANSFER BETWEEN TWO ID 
SEMICONDUCTORS 

The framework described in Sec. IHII can be easily 
generalized to electron transfer between two parallel ID 
semiconductors (intertube or interchain transfer). Real- 
izing its likely narrow region of applicability in terms of 
the offset of equilibrium energies (limited by the restric- 
tion to the lowest-energy electron states only), we will 
provide here just a brief outline. In what follows index i 
assumes values of 1 and 2 for the "donor" and "acceptor" 
systems. As the state of polarization is now described by 
two fictitious charge densities: Zi = J dxp{(x), the en- 
ergy stored in polarization is given now by 



U = 



\Y1 J dx d yp{( x ) v v( x -y)Pj(y)- 



Purely electronic states are derived from two Schrodinger 
equations 



Hi{x)^i{x) = (Ei-E^x) 
with their respective Hamiltonians: 

h 2 d 2 



(35) 



Hi{x) 



2m, dx 2 



%{x). 



Each is composed of the kinetic energy with the corre- 
sponding effective mass mi and the potential energy due 
to the solvent polarization: 



3 



dy V i:j (x-y) pUy) 



The lowest-energy solutions to Eqs. (|35j) yield the real 
charge densities Pi(x) = q\ipi(x)\ 2 . 

Following the same example of tubular charge distri- 
butions (radii Ri), the electrostatic kernels would be de- 
scribed by the following Fourier-transforms: 

V u {k) = 2Io(kR l )K (kR l )/e* 
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FIG. 4: Example of calculated electronic charge densities p(x) 
in the equilibrated polaron (solid line) and at the electron 
transition point (dashed line) for transfer between two iden- 
tical tubes. Parameters used here are R = 8 A, d = 19.5 A, 
and as = 32 A. 

for the intratube interactions (the same as Eq. (f!?Cj)) ). and 
V l2 {k) = 2I Q (kR 1 )I (kR 2 )K (kd)/e* (36) 

for the intertube interaction (instead of Eq. (|3"2")l ). where 
d is the distance between axes of the tubes. 

As in the earlier discussion, the unconstrained (in elec- 
tronic energy) minimization of total energies with the 
electron in state 1 or in state 2 results in the equilibrated 
polarons discussed in Sec. [IT] with some binding energies. 
The optimal electron transition state is determined under 
the constraint of equal electronic energies by minimizing 
the auxiliary functional 

£ = U + E l +a{E 2 -E l ) + y £ d (3 i Ni, (37) 

i 

where A/i is the normalization condition (|17p for wave- 
function tpi(x). One readily finds that the minimum of 
Eq. ([37]) is achieved for the fictitious densities satisfying 

p((x) = (1 - a)pi(x), p{{x) = ap 2 (x), 

that is, again obeying the partitioning (|25[) of the charge 
q at the transition point. One convenient approach to 
solving the problem of optimal fluctuations described by 
above equations is to specify the partitioning parameter 
a and then solve equations self-consistently to find the 
required energy offset and activation energies. 

A special interesting case is the transfer between two 
identical ID species, to which we restrict our comments 
here. In this case the optimal transition fluctuation cor- 
responds to a = 1/2. As an example of calculational 
results, Fig. [4] compares electronic density in the equili- 
brated polaron and at the electron transition point be- 
tween two identical tubes for a set of representative pa- 
rameters. As one should expect on physical grounds, the 



electronic density spreads out at the transition. Once 
again, the magnitude of the effect here is "mitigated" by 
the fact of relatively appreciable intertube electrostatic 
interaction (f3"6")) . From our calculations for a region of 
"physically reasonable" parameter values, we find that 
the activation energy in these cases is ordinarily a frac- 
tion (roughly 0.25-0.35) of the individual polaron binding 
energy Ef,. If the latter is on the order of 0.1 eV, this as- 
sessment would indicate that activation energies are com- 
parable to ambient temperatures making the thermally 
activated transfer between identical ID species efficient. 



V. DISCUSSION 

When a small-diameter ID semiconductor nanostruc- 
ture is immersed in a sluggish polar medium (solvent), 
fluctuations of the medium polarization may result in the 
appearance of localized electronic levels inside the band 
gap. An excess charge carrier (an electron or a hole) 
can occupy such a level and undergo self-localization into 
a large-radius adiabatic polaron accompanied by a self- 
consistent polarization pattern. Estimates based on rep- 
resentative parameter values indicate that the binding 
energy Ef, of solvation-induced ID polarons is expected 
to be on the order of 0.1 eV. Such estimates appear in 
agreement with our assessment [T3, El that the polaron 
binding may reach a substantial fraction, roughly one- 
third, of the binding energy of well-known Wannicr-Mott 
excitons in these structures, and experimental reports 
[H, MUl of the latter in SWCNTs as 0.4-0.6 eV. The 
physical properties of strong-coupling polarons are quite 
different from the nearly free band electrons frequently 
discussed in the context of electronic transport in nan- 
otubes and nanowires. One particularly notes that the 
mobility of solvated charge carriers is drastically reduced 
due to the dissipative drag of the medium [3, [l2|, l35| . 
On the other hand, thermal dissociation of excitons into 
electron- and hole-polaron pairs may be enhanced [10j . 

In this paper we addressed effects of the solvent po- 
larization on the fundamental process of electron trans- 
fer to/from ID semiconductors within an appropriately 
adapted framework of the Marcus theory of thermally 
activated transfer. The resulting qualitative picture 
shows that, for a range of energy offsets, the transfer 
can indeed take place via solvation-induced polaronic- 
like states whose energies are below the band continuum, 
at variance with the picture of transfer into rigid-band 
states unaffected by the polarization of the medium. In 
this sense, the transfer we discussed is closer to the tradi- 
tional scenario of transfer between two discrete electronic 
levels, both modulated by medium fluctuations. A qual- 
itative difference with the fixed-shape electronic states 
of small species, however, is that the electronic density 
on ID species is variable and can adjust its spatial ex- 
tent during the transfer. The framework we employed 
here considers only lowest-energy electronic states, and 
it would have to be extended with account of a multi- 
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tude of states in order to treat wider ranges of energy 
offsets and to verify if the inverted Marcus regime may 
be realizable. 

Similarly to approach of Ref. 25], our analysis was 
based on looking for optimal fluctuations of the medium 
polarization that promote iso-energetic electron transfer. 
Knowing the (free) energy of such fluctuations allows one 
to assess the corresponding thermal activation energies 
i? act . The optimal fluctuations correspond to the maxi- 
mum of the integrand in Eq. (J3J) and determine the expo- 
nential factor exp (—E^/kgT) in the transfer rate, ordi- 
narily of most interest in electron transfer problems [5j. 



The calculation of prefactors is more involved requiring 
considerations of polarization patterns "in the vicinity" 
of the optimal state at the transition. We refer the reader 
to Ref. [25[ for a discussion of some aspects of such cal- 
culations and just note here that the prefactors do not 
bring any "extraordinary" contributions 36]. 

Illustrative calculations in this paper have been per- 
formed for tubular ID geometries with possible applica- 
tions to carbon nanotube systems. The generic character 
of the picture should however be evident as applied to 
other transverse charge distribution geometries, such as 
in small-diameter nanowires or conjugated polymers. 
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